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The Schwinger-Dyson equations provide a convenient way to study dynamical 
symmetry breaking and dynamical mass generation in quantum field theory However, 
they constitute an infinite set of integral equations and surely some truncation scheme 
is necessary in order to be able to solve it. 

In pioneering works [|, the approach to dynamical fermion mass generation 
in quantum electrodynamics based on some truncated version of the Schwinger-Dyson 
equations (sometimes called ladder approximation) was developed, and it was shown 
the possibility of chiral symmetry breaking in QED, and also the existence of a critical 
coupling constant in the Landau gauge (for a review see ||). During the last fifteen 
years there has been much activity in developing and further extending the approach of 
refs. |], 3 to dynamical symmetry breaking in QED (for a review of the current status 
of this subject, see the proceedings 0). However, if dynamical mass is generated, the 
nonperturbative Ward-Takahashi identity in QED is not satisfied. That is why the 
dynamically generated fermionic mass, as well as the critical coupling constant in QED, 
are highly gauge-dependent. That fact was confirmed recently by the study of dynamical 
symmetry breaking in QED in an arbitrary covariant gauge ||. 

Currently, there is no doubt that the Schwinger-Dyson equations are powerful 
tools for the study of non-perturbative effects in field theory. However, perhaps some 
new approaches to truncation of the Schwinger-Dyson equations should be developed, 
particularly in order to solve the gauge- dependence problem in the dynamical symmetry 
breaking of QED. It may also be useful for this purpose to study the standard Schwinger- 
Dyson equations [fll in other models (albeit more complicated) like quantum gravity. 
Some time ago, a study along this line was done in Ref. || for AD Einstein quantum 
gravity coupled to fermions on a flat background. As a result, the possibility of chiral 
symmetry breaking in such model was shown ||. 

The purpose of the present letter is to study the Schwinger-Dyson equations and 
the problem of chiral symmetry breaking in 2D induced gravity coupled to fermions 
on a flat background. Notice that, recently, 2D induced gravity has been a popular 
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subject of study as a toy model for more realistic AD gravity, which still does not exist 
as a self-consistent theory. 

We shall start from the following action 



S - Sg + Sf, 



s a 



where R is the two-dimensional curvature, \l/ the 2D spinor, A the cosmological constant, 
and D^ is the 2D covariant derivative for spinors. In the standard approach to 2D 
induced gravity, working in the conformal gauge || [| 



(2) 



one can start from a matter theory (like Sf) and induce S g from it, by integrating over 
matter. After that, one has an exactly solvable theory with the action S g . 

Instead, we will adopt here another approach, namely, to start from the action (H), 
where gravitational and spinor fields are supposed to be quantum. Using the conformal 
gauge (0) and working in the vierbein formalism on a flat background, we consider 7 
as some given constant and do not integrate over Then, the term of the action 
corresponding to the interaction between \1/ and ip is given by 



s 



i nt 



- / d 2 x * 



(3) 



Let us now calculate the effective potential for composite fields [|T(J in the ladder 
approximation [0-0: 



V eS = -% Sp[lnVS - S^S + 1} + V 2 



(4) 



where the free fermion propagator is given by 



S (p) = -7, 



(5) 
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and the exact fermion propagator is 

s{p) = AW-Btfy (6) 

A and B being, for the moment, unknown functions. 

Note that V 2 corresponds (as it follows from the structure of S int (|3|)) to a two- 
particle irreducible diagram, which is similar to the analogous diagram in QED. Here, 
in V2 the fermion propagator is the exact one, while the vertex and the (p (graviton) 
field propagator 

are tree-level quantities (ladder approximation). Hence, V 2 is given by 

Using (|3|) and for calculating the effective potential (Ij), one can get, after per- 

forming Wick's rotation and the angular integration (we drop the details of these cal- 
culations): 



V, 







[fo dX 


In ^ 







B 2 (x)\ A(x)(A(x) - l)x + B 2 (x) 



x J xA 2 (x) + B 2 (x) 

+9 1" Tmhw) r ^whm [A{x)A(y)K4x ' y) + B ^ K ^ »>i} ■ 

(9) 

where Nf is the dimension of the fermion representation, M is the momentum cutoff, 
x = p 2 /M 2 , y = q 2 /M 2 , A(x) = A(p 2 ), B(x) = B(p 2 )/M and 



4xy + (x + y)(x + y + I - sj{x + y + I) 2 - Axy) 
K A (x,y) = 



2J(x + y + l) 2 - Axy 



2(x + y)+l- J(x + y + l) 2 - Axy 

K B (x,y) = -± y} Vl n / y - (10) 

(x + y + I) 2 — Axy 



A 



with the notations 
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" 64:71' 

and 



'11' 



1 = 2^' < 12 > 
The Schwinger-Dyson equation which corresponds to the effective potential of 

the composite fields (|9]), in the ladder approximation takes the form 

72 

l (S- 1 (p)-S 1 (p))= J — ^ [r(p - q, q)S{q)T{q - p,p)G(p, q)} . (13) 

Bearing in mind integrating over the angles and doing the same changes of variable 
as before in ( |T3"D , one can show that the functions A and B must obey integral equations 
of the following form 

A(x) = l+ 'Sl \^mST ' KA{t ' v) - 

B(x) = ' £ \A>w + B>M K, * x 'ti- (14) 

Of course, it is not possible to solve these equations analytically. However, we will 
obtain their numerical solution by a standard iterative method at some region of the 
theory parameters. First, we fix the values of the parameters g and I appearing in (|T2|) 
and prepare two types of trial functions (the procedure is quite similar to the one used 
inref. @): 

(a) A°(x) = ci, B°(x) = 0, 

(b) A°(x) = Cl , B°(x) = c 2 , 

where C\ abd c 2 are constants between and 1. The functions A°(x) and B°(x) are 
the starting point of a self-consistent iterative calculation in which one finds successive 
pairs of functions A l (x) and B l (x). To be more precise, each pair is obtained from the 
previous one by means of the recurrence 

r i A i (v) 1 

s ' +1(i) = »I?* wS^) Kb(x ' v) - (15) 



The sequences formed by the {A* (a;)} and {B*(x)} are expected to converge into the 
functions A(x) and B(x), respectively, which are the desired solutions of fll4|) . In prac- 
tice, we judge the degree of convergence of these series by the smallness of the squared 
norms of the differences A l+1 — A 1 and B l+l —B\ In our calculations, we have set bounds 
of orders 10~ 4 — 10~ 6 . If for the given g and / there are solutions of both types, (a) and 
(b), only the most stable of both by V e s © must be chosen as the one corresponding 
to the true vacuum. 

We have executed this algorithm to solve (0), starting from trial functions (a) 
and (b) for fixed I = 0.5 and varying g. For very small g's, both types lead to curves 
close to A(x) = 1, B(x) = 0, i.e. the chiral symmetric solution, as should be expected. 
Moreover, their respective potentials are practically undistinguishable. As g increases, 
the Kff for the (a)-type solution appears to be just slightly higher than for the (b)-type 
one, thus selecting the latter as the physical vacuum. This happens until some specific 
value of g, around 0.1, is reached, for which the (a)-solution looks rather reluctant to 
converge. Before this happens, V e s is marginally higher than for the (b)-solution. This 
seems to indicate that the solution 'lost' was not physical. 

Typical curves representing the A and B functions obtained are shown in Fig. 1. 
All of them correspond to (b)-solutions (i.e. chiral non-symmetric solutions) found for 
different values of g, and for I = 0.5. The starting constants where c\ = C2 = 0.5, but 
their precise value has no influence on the final form of A and B. 

As is clear, the value of V e g corresponding to a solution obtained for a given g 
changes as we vary g. V e s vanishes as g approaches zero, and increases as g grows. 
There is a specific g for which it stops growing, and shows a slow decrease for larger g's. 
However, this should be regarded cautiously, since the larger g is, the closer we are to 
regions where perturbative methods fail. 

Hence, as we see from our numerical analysis, there is a possibility of chiral 
symmetry breaking in 2D induced gravity coupled to fermions. Our study has been 
done in the 'physical' conformal gauge which, in some sense, may be considered as 
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the analogue of the 'physical' Landau gauge in QED. The results of these numerical 
studies should be gauge dependent, particularly because we do not pay attention to the 
nonperturbative Ward-Takahashi identity resulting from general covariance. (Note that 
this identity is not so easy to obtain explicitly in the noncovariant conformal gauge). 
Hence, for a complete proof that chiral symmetry breaking in 2D induced gravity is a 
real physical phenomenon one has to repeat the same study in some covariant gauge (of 
harmonic type 11], for example). This is left to be done elsewhere. 
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Figure Caption 

Fig. 1. 

Plot of the functions A and B obtained as the (b)-type solutions for g =0.1, 0.2 
and 0.25, keeping / = 0.5 fixed. Notice how B deviates more and more from the g = 
solution (B(x) = 0) as g increases. Although not shown in the figure, the curve keeps 
going up for larger values of g. 
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